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ABSTRACT 

Cosmological implications of the observed large-scale peculiar velocities are 
reviewed, alone or combined with redshift surveys and CMB data. The latest ver- 
sion of the POTENT method for reconstructing the underlying three-dimensional 
velocity and mass-density fields is described. The initial fluctuations and the na- 
ture of the dark matter are addressed via statistics such as bulk flow and mass 
power spectrum. The focus is on constraining the mass density parameter Q, 
directly or via the parameter f3 which involves the unknown relation between 
galaxies and mass. The acceptable range for Q is found to be 0.4 — 1.0. The 
range of jS estimates is likely to reflect non-trivial features in the galaxy biasing 
scheme, such as scale dependence. Similar constraints on Q and A from global 
measures are summarized. 



1. INTRODUCTION 

A major goal of the analysis of cosmic flows is measuring the cosmological parameters, 
in particular the mass density parameter, Q (or Q). As illustrated in Figure |l], the data can 
be used to constrain Q in several different ways. Methods that are based on the peculiar 
velocity data alone (§ ^ are independent of the "biasing" relation between galaxy and mass 
density. They refer to the present-day large-scale structure (LSS) which is insensitive to 
the cosmological constant A. They can thus serve to measure Q directly, but they involve 
relatively large errors. These methods largely rely on the assumption (supported by obser- 
vations, e.g., Nusser, Dekel & Yahil 1995) that the initial fluctuations were drawn from a 
Gaussian random field. 

All the methods that use the spatial distribution of galaxies must depend on the biasing 
relation between the densities of galaxies and mass. In the linear approximation to gravita- 
tional instability theory (GI), which is roughly valid for the fields when they are smoothed 
on very large scales, what is actually being measured by these methods is the degenerate 



^in "Galaxy Scaling Relations: Origins, Evolution and Applications", ed. L. da Costa (Springer) in press 
(1997). A more complete discussion is in "Formation of Structure in the Universe", eds. A. Dekel & J. P. 
Ostriker (Cambridge University Press) in press (1997). 
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Fig. 1. — Methods for measuring 17 and (3 from large-scale flows. 

parameter (3 = Q^'^/b rather than Q itself, where the density fluctuations of galaxies and 
mass are assumed to be related via a linear biasing relation, 6g = b6. These methods in- 
clude, for example, measurements of redshift-space distortions from redshift surveys under 
the assumption of global isotropy, comparisons of the galaxy distribution with the CMB 
dipole, and comparisons of the observed peculiar velocities with the galaxy distribution (or 
the predicted velocities) deduced from redshift surveys (§ The most recent best estimates 
of f3 for IRAS galaxies lie in the range 0.5 < (3i < 0.9. As argued below, the contamination 
by non-trivial biasing introduces a significant uncertainty in the translation of the various j3 
estimates to a reliable measurement of Vt. 

Another way of determining VL is by comparing the mass power spectrum from observed 
peculiar velocities with the spectrum of fluctuations in the CMB (§ pl3| , § p.4|) . This compari- 
son involves more detailed modeling of the formation of structure, and the current studies are 
limited to the CDM family of models, allowing for a non-zero cosmological constant, a pos- 
sible tilt in the spectrum with or without tensor fluctuations, and a mixture of dark-matter 
species. 

The outline is as follows: § ^ briefly describes reconstruction methods from peculiar 
velocities. § ^ discusses the statistics of mass-density fluctuations, including bulk flow and 
power spectrum, as determined by the velocity data alone or combined with CMB data. § ^ 
focuses on direct estimates of VL from peculiar velocities alone. § ^ presents estimates of (3 
by comparing velocity and galaxy density data, and addresses the issue of biasing. § ^ puts 
the results in perspective by summarizing other measures of the cosmological parameters. 
§ |6.1| provides a summary of the results. 



2. POTENT RECONSTRUCTION FROM PECULIAR VELOCITIES 
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2.1. Data for Velocity Analysis 



The key for measuring peculiar velocities are the distance indicators (for review: Willick 
1997). The current typical intrinsic scattering the TuUy-Fisher (TF) distance indicator is at 
best (Jm~0.33 mag, corresponding to a relative distance error of A= (In 10/5)crm~0.15. 

The most comprehensive catalog of peculiar velocity data available today is the Mark III 
catalog (Willick et al. 1995; 1996; 1997a), which is a careful compilation of several data sets 
under the assumption that all galaxies trace the same underlying velocity field. The merger 
was non-trivial because the observers differ in their selection procedure, the quantities they 
measure, the method of measurement and the TF calibration techniques. 

The original Mark II catalog, which was used in the first application of POTENT (Dekel 
et al. 1990; Bertschinger et al. 1990), consisted of about 1000 galaxies (mostly Lynden-Bell 
et al. 1988 and Aaronson et al. 1982). The extended Mark III catalog consists of ~ 3400 
galaxies (dominated by Mathewson et al. 1992). This sample enables a reasonable recovery 
of the dynamical fields with 12/;,~^Mpc smoothing in a sphere of radius ~60/z~^Mpc about 
the Local Group (LG), extending to ~80 /i~^Mpc in certain regions (§ |2.7|) . More uniformly 
sampled data of more than 1000 spiral galaxies in the north is in preparation, and was 
already subject to preliminary analysis (SFI, Giovanelli et al. 1997). 



2.2. Methods of Velocity Analysis 

One way of classifying the methods of velocity analysis is as follows: 





Inferred Distance Space 


Redshift-Space + V Model 


Forward 
TF 


POTENT 

selection + distance bias 


VELMOD 

selection bias 


Inverse 
TF 


POTINV 

distance bias 


MFPOT 

z-space smoothing 



The "forward" and "inverse" methods refer to whether the TF relation is interpreted as M{ri) 
or ri{M) (M being the absolute magnitude and t] the rotation velocity). The difference is 
crucial because the apparent magnitude depends on distance while r] is not, and because the 
selection depends on magnitude and is independent of r]. On the other hand, the velocity 
field can be computed either from peculiar velocities that were evaluated at each galaxy's 
TF-inferred position, d, or by fitting a parametric model for the potential (and thus the 
velocity) field in redshift space, z. Each method is affected by different systematic errors, 
and techniques have been developed for statistically correcting them. The success of each 
technique has been tested using mock catalogs, and the goal is to have the different methods 
recover consistent results. 

The original POTENT described below is a forward TF method in d-space. It has to 
deal with selection bias of the TF parameters because of the magnitude limit, and with 
Malmquist bias in the inferred distances and velocities arising from the random distance 
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errors convolved with the geometry of space and the dumpiness of the galaxy distribution. 
Methods have been developed for statistically correcting these biases. In particular, the 
correction of the Malmquist bias requires external information about the underlying number 
density of galaxies in the samples from which galaxies were selected for the peculiar velocity 
catalogs (see Willick 1997). 

One can alternatively infer distances using the inverse TF relation. This eliminates the 
selection bias, but there is still an inferred-distance Malmquist bias. In the inverse- TF case, 
the distance bias can in principle be corrected using information that is fully contained in 
the catalog itself (Landy & Szalay 1992). In practice, the quality of the correction is limited 
by the sparseness of the sampling. The POTENT analysis of the inverse data corrected this 
way is termed POTINV (Eldar, Dekel & Willick 1997). 

If the selection does not explicitly depend on r], Malmquist bias can be eliminated by 
minimizing rj residuals in redshift space without ever inferring actual distances to individual 
galaxies. The distance is replaced hj r = z — Ua{z), where Ua is a parametric model for the 
radial peculiar velocity field. If the forward TF relation is used, as in VELMOD (Willick et 
al. 1997b; § |5.4|) , the method still has to correct for selection bias. The use of the inverse 
TF relation (Schechter 1980) guarantees in this case that both the selection bias and the 
distance bias are eliminated, at the expense of over-smoothing due to the representation of 
the fields in redshift space. Recent implementations of such methods were termed MFPOT 
(Blumenthal, Dekel & Yahil 1997) and ITF (Davis, Nusser, & Willick 1996). 

Another way of distinguishing between the methods is by their goals. The methods 
working in d-space can serve for reconstruction of 3D maps of the velocity and mass-density 
fields, unbiased and uniformly smoothed with equal-volume weighting throughout the vol- 
ume. These fields can then be straightforwardly compared to other data and to theory in 
order to obtain cosmological implications {e.g., § [5.2| ). Alternatively, one may direct the 
method to estimating certain parameters of the model {e.g., (3) without ever reconstructing 
uniform maps. The redshift-space methods serve this purpose well {e.g., § [5.4| ). 

Yet another characteristic of some of the methods is the usage of a whole-sky redshift 
survey (such as IRAS 1.2 Jy) as an intrinsic part of the reconstruction from peculiar velocities. 
This is the case in the SIMPOT, VELMOD and ITF methods (§ These methods are 
geared towards determining (3, with SIMPOT also providing uniform reconstruction maps. 

Finally, one can focus on optimal formal treatment of the random errors, which are 
in fact the main obstacle. A method based on Wiener Filtering has been developed for 
recovering the most probable mean field from the noisy peculiar-velocity data in d-space 
(Zaroubi, Dekel & Hoffman, in preparation). This can serve as a basis for constrained 
realizations of uniform smoothing, each of which being an equally good guess for the structure 
in our real cosmological neighborhood. 
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2.3. Correcting Malmquist Bias 

The selection bias in the calibration of the forward TF relation can be corrected once 
the selection function is known (see Willick 1994). But then, the TF inferred distance, and 
the mean peculiar velocity at a given d, suffer from a Malmquist or inferred- distance bias. 
The distances, either forward or inverse, are corrected for Malmquist bias in a statistical way 
before being fed as input to POTENT-like procedures. 

If M is distributed normally for a given 77, with standard deviation 0"^, then the forward 
inferred distance d oi a. galaxy at a true distance r is distributed log-normally about r, with 
relative error A^0.46(T.m- Given d, the expectation value of r is 



S^P{r\d)dr J^^rMr) exp {~^^) dr 



where n(r) is the number density in the underlying distribution from which galaxies were 
selected. The deviation of E{r\d) from d reflects the bias. The homogeneous part arises 
from the geometry of space — the inferred distance d underestimates r because it is more 
likely to have been scattered by errors from r > d than from r < d, the volume being oc r^. 
If n = const, equation ^ reduces to E{r\d) = de^'^^ , in which the distances should simply 
be multiplied by a factor, 8% for A = 0.15, equivalent to changing the zero-point of the 
TF relation. Fluctuations in n(r) are responsible for the inhomogeneous bias (IM), which 
systematically enhances the inferred density perturbations and thus the value of Q inferred 
from them. 

In one version of the Mark III data for POTENT analysis, the forward IM bias is 
corrected in two steps. First, the galaxies are grouped in z-space (Willick et al. 1995), 
reducing the distance error of each group of members to A/^/N and thus significantly 
weakening the bias. With or without grouping, the noisy inferred distance of each object, 
d, is replaced by E{r\d), with an assumed n{r) properly corrected for grouping if necessary. 
This procedure has been tested using realistic mock data from N-body simulations, showing 
that IM bias can be reduced to the level of a few percent. The practical uncertainty is in 
n{r), which can be approximated for example by the high- resolution density field of IRAS 
or optical galaxies, or by the recovered mass-density itself in an iterative procedure under 
certain assumptions about how galaxies trace mass. The resultant correction to the density 
recovered by POTENT is < 20% even at the highest peaks. 

Distances are alternatively inferred via the inverse TF relation between internal velocity 
parameter t] and magnitude m, t] = ri^{m — 5\ogd). Under the assumption that the selection 
was independent of f] and was not an explicit function of distance, the expectation value of 
the true distance r given d is 

E{r\d) = dc'^'/'f{dc^')/f{d) , (2) 

where A = (In 10/5)cr^/r7°. In this case, the required density function, f{d), is in d-space, 
and is derivable from the sample itself (Landy & Szalay 1992). Eldar, Dekel & Willick 
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(1997) have applied this correction to the inverse distances in the Mark III catalog, to serve 
as input for a POTENT analysis (POTINV). The agreement between the forward POTENT 
and POTINV results are well within the level of the random errors. 



2.4. Smoothing the Radial Velocities 

The goal of the POTENT analysis is to recover from the collection of Malmquist- 
corrected, radial peculiar velocities Ui at inferred positions di the underlying 3D velocity 
field v{x) and the associated mass-density fluctuation field S{x), smoothed with a Gaussian 
of radius Rg (we denote hereafter a 3D Gaussian window of radius 12 h~^Mpc by G12, etc.). 
The first, most difficult step is the smoothing, or interpolation, into a radial velocity field 
with minimum bias, u{x). The desire is to reproduce the u{x) that would have been ob- 
tained had the true v{x) been sampled densely and uniformly and smoothed with a spherical 
Gaussian window of radius Rs- With the data as available, u{Xc) is taken to be the value 
at x = Xc of an appropriate /oca/ velocity model v{ak^x — Xc). The model parameters are 
obtained by minimizing the weighted sum of residuals, 

S = Y,Wi[ui-Xi-v{ak,Xi)f , (3) 

i 

within an appropriate local window Wi = W{Xi,Xc)- The window is a Gaussian, modified 
such that it minimizes the combined effect of the following three types of errors. 

Tensor window bias. Unless Rg r, the UjS cannot be averaged as scalars because the 
directions Xi differ from Xc, so u{Xc) requires a fit of a local 3D model as in Eq. ^ The original 
POTENT used the simplest local model, v(x) =5 of 3 parameters, for which the solution can 
be expressed explicitly in terms of a tensor window function (Dekel et al. 1990). However, a 
bias occurs because the tensorial correction to the spherical window has conical symmetry, 
weighting more heavily objects of large Xi-Xc- A way to reduce this bias is by generalizing 
the zeroth-order B into a 9-parameter first-order velocity model, v{x) =B+L-{x—Xc), with 
L a symmetric tensor that automatically ensures local irrotationality. The linear terms tend 
to "absorb" most of the bias, leaving v{Xc) = B less biased. Unfortunately, a high-order 
model tends to pick undesired small-scale noise. The optimal compromise for the Mark III 
data was found to be a 9-parameter model fit out to r = 40 /i~^Mpc, smoothly changing to a 
3-parameter fit beyond 60 /i~^Mpc (Dekel et al. 1997). 

Sampling-gradient bias. If the true velocity field is varying within the effective win- 
dow, the non-uniform sampling introduces a bias because the smoothing is galaxy-weighted 
whereas the aim is equal-volume weighting. The simplest way to correct this bias is by 
weighting each object with the local volume that it "occupies", or the inverse of the local 
density. A crude estimate of this volume is oc i?^, where i?„ is the distance to the n- 
th neighboring object {e.g., n = 4). This procedure is found via simulations to reduce the 
sampling-gradient bias in Mark III to negligible levels typically out to 60 /i^^Mpc as long as 
one keeps out of the Galactic zone of avoidance. The Rn{x) field can serve later as a fiag for 
poorly sampled regions, to be excluded from any quantitative analysis. 
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Reducing random errors. The ideal weighting for reducing the effect of Gaussian noise 
has weights Wi(xa~'^, where are the distance errors. Unfortunately, this weighting spoils 
the carefully designed volume weighting, biasing u towards its values at smaller r, and at 
nearby clusters where the errors are small. A successful compromise is to weight by both, 
i.e. 

W{Xi,Xc) oc Vi cr-' e^p[-{Xi - Xcf/2Rl] . (4) 

The resultant errors in the recovered fields are assessed by Monte-Carlo simulations. 
We generate noisy data via full, realistic Monte-Carlo mock catalogs, where the noise is 
added as scatter in the TF quantities (Kolatt et al. 1996). The error in the final 5 at a 
grid point is estimated by the standard deviation of the recovered 5 over the Monte-Carlo 
simulations, (and similarly cr„). In the well-sampled regions, which extend in Mark III 
out to 40 — 60/i~^Mpc, the errors are cr^Ri 0.1 — 0.3, but they may blow up in certain regions 
at large distances. To exclude noisy regions, any quantitative analysis could be limited to 
points where Uy and cr^ are within certain bounds. 



2.5. From Radial Velocity to Density Fields 

If the LSS evolved according to Gl, then the large-scale velocity field is expected to be 
irrotational, Vxv — 0. Any vorticity mode would have decayed during the linear regime as the 
universe expanded, and, based on Kelvin's circulation theorem, the fiow remains vorticity- 
free in the mildly-nonlinear regime as long as it is laminar. Bertschinger & Dekel (1989) have 
demonstrated that irrotationality is valid to a good approximation when a nonlinear velocity 
field is properly smoothed over. Irrotationality implies that the velocity field can be derived 
from a scalar potential, v{x) — — V$(a:), so the radial velocity field u{x) should contain in 
principle enough information for a full 3D reconstruction. In the POTENT procedure, the 
potential is computed by integration along radial rays from the observer, 

^x) = - r u(r',e,(j))dr' . (5) 
Jo 

The two missing transverse velocity components are then recovered by differentiation. 

The final step of the POTENT procedure is the derivation of the mass-density fiuctua- 
tion field associated with the peculiar velocity field. This requires a solution to the equations 
of GI in the mildly-nonlinear regime with mixed boundary conditions. 

Let X, V be the position and peculiar velocity in comoving units (corresponding to ax 
and av in physical units, with a{t) the universal expansion factor). Let 5 = {p — p)/p he 
the mass-density fiuctuation. The equations governing the evolution of fiuctuations of a 
pressureless gravitating ffuid in a standard cosmological background during the matter era 
are the Continuity equation, the Euler equation of motion, and the Poisson field equation. 

In the linear approximation, the growing mode of the solution, 5 oc D(t), is irrotational 
and can be expressed in terms of f{Q,) = H~^D/D^QP-^. The corresponding hnear relation 
between density and velocity is 5i — —f~^V-v. The use of 5-i is hmited to the small 
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dynamical range between a few tens of megaparsecs and the ~ 100 h~^Mpc extent of the 
current samples. However, the sampling of galaxies enables reliable dynamical analysis 
with a smoothing radius as small as ~ 10 /i~^Mpc, where |V-w| obtains values larger than 
unity and therefore nonlinear effects play a role. Even reconstruction with ~ 5 h~^Mpc 
smoothing may be feasible in well-sampled regions nearby. Figure |^ shows that 6i becomes 
a severe underestimate at large \5\. Mild nonlinear effects carry crucial information about 
the formation of LSS, and should therefore be treated properly. 
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Fig. 2. — Quasi-linear velocity-to-density approximations. A5 = Sapprox{v) — Stme- The mean 
and standard deviation are from large CDM N-body simulations normalized to fig = 1, smoothed 
G12. The three curves correspond to different models: standard CDM (squares, solid), tilted CDM 
(n = 0.6, stars, dashed), and open CDM {Q = 0.2, open circles, dotted) (Canon et al. 1997). 

A basis for useful mildly-nonlinear Telations is provided by the ZeVdovich (1970) approx- 
imation. The displacements of particles from their initial, Lagrangian positions q to their 
Eulerian positions x at time t are assumed to have a universal time dependence, and thus, 
^(9) ~ 9 = For the purpose of approximating GI, the Lagrangian Zel'dovich 

approximation can be interpreted in Eulerian space, q{x) =x — f~^v{x), provided that the 
flow is laminar with no orbit mixing, or when multi-streams are appropriately smoothed 
over. The solution of the continuity equation then yields (Nusser et al. 1991) 

6c{x) = \\I - f-^dv/dx\\-l , (6) 

where the bars denote the Jacobian determinant and / is the unit matrix. The Zel'dovich dis- 
placement is flrst order in and v and therefore the determinant in 6c includes second- and 
third-order terms as well, involving sums of double and triple products of partial derivatives 
which we term A2(a;) and A^lx) respectively. 
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The approximation 5c can be improved by slight adjustments of the coefficients of the 
nth-order terms, 

5,+ = -(1 + tr)f-^V-v + (1 + e2)/-'A2 + (1 + e3)/-^A3 . (7) 

The coefficients were empirically tuned to best ffi the CDM simulation of 12 /i^^Mpc smooth- 
ing over the whole range of 5 values, with ei = 0.06, e2 = —0.13 and €3 = —0.3. This approx- 
imation is found to be robust to uncertain features such as the value of fi, the shape of the 
power spectrum, and the degree of non-linearity as determined by the fluctuation amplitude 
and the smoothing scale. Such robustness is crucial when a quasilinear approximation is 
used for determining f2 (§ H). This is the approximation currently used in POTENT. 

Fig. 1^ compares the accuracy of the various approximations using the N-body simula- 
tions. 5c is the best among the physically motivated approximations, which also include two 
second-order approximations (Bernardeau 1992; Gramman 1993). The latter do somewhat 
better at the negative tail, but they provide severe underestimates in the positive tail. 5c+ 
is an excellent robust fits over the whole mildly-nonlinear regime. 

We note in passing that the relation ^ is not easily invertible to solve for v when 5 is 
given, e.g., from redshift surveys, but a useful approximation derived from simulations is 
V-t; = + 0.18(5). 



2.6. Testing with Mock Catalogs 

The way to optimize POTENT and other reconstruction methods is by minimizing the 
systematic errors when applied to mock catalogs. It is important that these mock catalogs 
mimic the real data as closely as possible. Such mock catalogs have been produced, for 
example, to mimic the Mark III and the IRAS 1.2Jy catalogs (Kolatt et al. 1996). They are 
publically available and serve as standard "benchmarks" for the competing methods. 

The procedure for making these mock catalogs involves two main steps: a dynamical 
N-body simulation that mimics our actual cosmological neighborhood, and the generation 
of galaxy catalogs from it. 

Figure |^ demonstrate the quality of the POTENT reconstruction from the Mark III 
catalog by comparing the recovered density field to the true G12-smoothed field. This com- 
parison is done at the points of a uniform grid inside a volume of effective radius 40 /;,~^Mpc. 
The field shown is the average of the fields recovered from ten Monte Carlo mock catalogs of 
noisy galaxy velocities sampled sparsely and nonuniformly. One can see that the remaining 
systematic errors are small. The final systematic error is not correlated with the signal (slope 
~unity in the scatter diagram) and is on the order of A5 ~ 0.13. The random errors are 
not a major obstacle in certain well-sampled regions (such as the Great Attractor), but they 
become severe in poorly-sampled regions (such as parts of the Perseus-Pisces region near 
the Galactic plane). The errors derived from the noisy mock catalogs are used to eliminate 
poorly-recovered regions from quantitative analyses. 
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Fig. 3. — Systematic errors in tlie POTENT analysis. The density field recovered by POTENT 
from the mock data is compared with the "true" G12 density. The comparison is at uniform 
grid points within a volume of radius 40/i~^Mpc. Left: The input to POTENT is the true, G12- 
smoothed radial velocity. The small scatter of 2.5% reflects small deviations from potential flow, 
the scatter in the non- linear approximation Eq. and numerical errors. Right: The input is noisy 
and sparsely-sampled mock data. Shown is the average of 10 random realizations. The global bias 
is only -4% (Dekel et al. 1997). 

2.7. Maps of Velocity and Density Fields 

Figure ^ shows Supergalactic-plane maps of the velocity field in the CMB frame and 
the associated 5c+ field (for f2 = 1) as recovered by POTENT from the Mark III catalog. The 
recovery is reliable out to ~60 /i~^Mpc in most directions outside the Galactic plane (F = 0). 
Both large-scale (~ 100 /i^^Mpc) and small-scale (~10/i~^Mpc) features are important; e.g., 
the bulk velocity reflects properties of the initial fluctuation power spectrum (§ |^), while the 
small-scale variations indicate the value of VL (§ |^). 

The velocity map shows a clear tendency for motion from right to left, in the general di- 
rection of the LG motion in the CMB frame (L, i? = 139°, —31° in Supergalactic coordinates). 
The bulk velocity within 60 /I'^Mpc is 300 - 350 km s"^ towards (L, 5^166°, -20°) (§ |3l|) 
but the flow is not coherent over the whole volume sampled, e.g., there are regions in front of 
PP (bottom right) and behind the GA (far left) where the XY velocity components vanish, 
i.e., the streaming relative to the LG is opposite to the bulk flow direction. The velocity 
field shows local convergences and divergences which indicate strong density variations on 
scales about twice as large as the smoothing scale. 

The bottom panel of Fig. ^ shows the POTENT mass density field in the Supergalactic 
plane as a landscape plot. The Great Attractor (with 12 /I'^Mpc smoothing and ^2 = 1) is a 
broad density ramp of maximum height 5 = 1.4 ±0.3 located near the Galactic plane y = at 
X f» — 40 /i~^Mpc. The GA extends towards Virgo near F ^ 10 (the "Local Supercluster"), 
towards Pavo- Indus- Telescopium (PIT) across the Galactic plane to the south (F <0), and 




Fig. 4. — POTENT dark-matter maps. The G12 fluctuation fields of peculiar velocity and mass- 
density in the Supergalactic plane as recovered by POTENT from the Mark III peculiar velocities. 
The vectors are projections of the 3D velocity field in the CMB frame, dominated by large deriva- 
tives embedded in a coherent bulk flow. Distances and velocities are in lOOkms"^. Contour spacing 
is 0.2 in 6, with the heavy contour marking 5 = and dashed contours 6 < 0. The height in the 
surface plot is proportional to the total mass density contrast 6. The LG is at the center, GA on 
the left, PP and the Southern Wall on the right. Coma Great Wall at the top, and the Sculptor 
void in between (Dekel et al. 1997). 

towards the Shapley concentration behind the GA (F > 0,X < 0). The structure at the 
top is related to the "Great Wall" of Coma, with 5 ~ 0.6. The Perseus Pisces peak which 
dominates the right-bottom is peaked near Perseus with 5 = 1.0 ± 0.4. PP extends towards 
the southern galactic hemisphere (Aquarius, Cetus), coinciding with the "Southern Wall" 
as seen in redshift surveys. Underdense regions separate the GA and PP, extending from 
bottom-left to top-right. The deepest region in the Supergalactic plane, with 5 = — 0.8 ±0.2, 
roughly coincides with the galaxy- void of Sculptor and is useful in bounding VL (§ §). 
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3. STATISTICS OF MASS-DENSITY FLUCTUATIONS 

Having assumed evolution by GI, the structure can be traced backward in time in order 
to recover the initial fluctuations and to measure statistics which characterize them as a 
random field, e.g., the power spectrum, P{k), and the probability distribution functions 
(PDF). "Initial" here may refer either to the linear regime at 2;~10^ after the onset of the 
self-gravitating matter era, or to the origin of fiuctuations in the early universe before being 
filtered on sub-horizon scales during the plasma-radiation era. The spectrum is filtered 
on scales < 100 h~^Mpc by dark-matter dominated processes, but its shape on scales > 
10 h~^Mpc is not affected much by the mildly-nonlinear effects (because the faster density 
growth in superclusters roughly balances the slower density depletion in voids at the same 
wavelength). The shape of the one-point PDF, on the other hand, is expected to survive the 
plasma era unchanged but it develops strong skewness even in the mildly-nonlinear regime. 
Thus, the present day P{k) can be used as is to constrain the origin of fluctuations (on large 
scales) and the nature of the DM (on small scales), while the PDF needs to be traced back 
to the linear regime flrst. 

The competing scenarios of LSS formation are reviewed for example by Primack (1997). 
In summary, if the dark matter (DM) is all baryonic, then by nucleosynthesis constraints 
the universe must be of low density, f2 ^ 0.1, and a viable model for LSS is the Primordial 
Isocurvature Baryonic model (PIB) with several free parameters. With f2 ~ 1, the non- 
baryonic DM constituents are either "cold" or "hot", and the main competing models are 
CDM, HDM, and CHDM — a 7:3 mixture of the two. The main difference in the DM 
effect on P{k) arises from free-streaming damping of the "hot" component of fluctuations 
on galactic scales. Currently popular variants of the standard CDM model (fi = 1, n = l) 
include a tilted power spectrum on large scales [n ^ 1) and a flat, low-f2 universe with a 
non-zero cosmological constant such that Q + fl\ = 1. 

The peculiar velocities of the Mark III catalog enable direct derivations of the mass 
power spectrum itself, independent of galaxy biasing, roughly in the range 10—100 h~^Mpc. 
The bulk velocity in spheres of radii up to 60 h~^Mpc is sensitive to even larger wave- 
lengths. In all standard theories, the power spectrum on large scales is expected to be 
a power law, Pk oc fc", with n of order unity. It is expected to turn around at kpeak ~ 
0.065(r2/;.)^^( /i^^Mpc)"^, corresponding to the horizon scale at the epoch of equal energy 
densities in matter and radiation. The dark matter type mostly affects the shape of the 
flltered spectrum in the "blue" side of the peak {k > kp^ak)- Once the fluctuation amplitude 
on very large-scales is flxed by COBE's measurements of CMB fluctuations, the bulk velocity 
is sensitive to n and is insensitive to fl or the DM type. The steep slope of the CDM-like 
spectra at A; > kp^ak, where it is best constrained by the data, makes it more sensitive than 
the bulk velocity to Qh. 

We flrst describe the bulk velocity (§ p.l|) . Then two ways of evaluating P{k): a model- 
independent evaluation from the velocity fleld via POTENT (§ p.2|) , and a likelihood esti- 
mation from raw radial peculiar velocities under a prior model (§ p. 31 ). The P{k) from the 
local velocities is then compared to sub-degree angular power spectrum of CMB fluctuations 
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(§0). 



3.1. Bulk Velocity 

A simple and robust statistic related to the power spectrum is the bulk velocity — the 
amplitude of the vector average V of the i?s-smoothed velocity field v{x) over a volume 
defined by a normalized window function WR{r) {e.g., top-hat) of a characteristic scale R, 

V = j d'x Wr{x) v{x) , ly'') = dk P{k) W^k) . (8) 

We denote by Vr the bulk velocity in a top-hat sphere of radius R = r h~^Mpc. The ensemble 
variance (V^^) for a model that is characterized by P{k) is an integral of P{k) in which the 
wavelengths > R are emphasized by W^{k), the Fourier transform of WR{r). The bulk 
velocity can be obtained from the observed radial velocities by minimizing residuals as in 
Eq. ^ The first report by Dressier et al. (1987) oi V = 599 ± 104 for ellipticals within 
~60/i~^Mpc was interpreted prematurely as being in severe excess of common predictions, 
but it quickly became clear that the effective window was much smaller due to the nonuniform 
sampling and weighting (Kaiser 1988). The sampling-gradient bias can be crudely corrected 
by volume weighting as in POTENT (§ |2.4| ), at the expense of larger noise. Courteau et 
al. (1993) reported based on an early version of the Mark 111 data Veo = 360 ± 40 towards 
(L,S) = (162°, -36°). Alternatively V, can be computed from the POTENT v field by 
simple vector averaging from the grid. 

The bulk velocity as a function of i?, from several recent sources, is shown in Figure 
|. The Mark III POTENT result ai R = 50 /I'^Mpc is V^o = 374 ± 85kms-^ towards 
(158°, —9°) ± 10°. The ~20% error bars are due to distance errors, and one should consider 
an additional uncertainty of similar magnitude due to the non-uniform sampling. The SFI 
sample of Sc galaxies yields at the same R a very similar result (contrary to premature 
rumors), V50 ~ 364kms~^ towards (172°, —14°) (daCosta et al. 1996). These samples are 
not large enough for a reliable estimate of V at larger radii. 

Supernovae Type la provide more accurate distances, with only ~ 8% error, and they 
can be measured at larger distances. The current sample of 44 such SNe by Riess & Kirshner 
(1997, following Riess, Press & Kirshner 1995), which extends out to ~ 300/;,~^Mpc, shows 
a bulk flow of V = 253 ± 252kms~^ towards (166°, —44°). The effective radius of this data 
set for a bulk flow flt is in fact less than 50 h^^Mpc because the data is weighted inversely 
by the errors. The SNe bulk flow is consistent with the results from the Mark III and SFI 
galaxy data. They all make a bulk of sense within the framework of standard isotropic and 
homogeneous cosmology. 

The only apparently discrepant result comes from the velocities measured on a larger 
scale using brightest cluster galaxies (BCG) as distance indicators (Lauer & Postman 1994, 
LP). They indicate a large bulk velocity of = 689 ± 178 towards a very different direction 
~ (126, 21). An ongoing effort to measure BCG's in a larger sample of clusters and distances 
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Fig. 5. — Bulk velocity. The amplitude of the bulk velocity relative to the CMB frame in top- 
hat spheres about the LG, as derived from several data sets. The directions are indicated in 
supergalactic coordinates {L,B). The Mark III data and the SFI data yield consistent results. The 
new result from Supernovae type la on larger scales is a natural extrapolation. The result from 
brightest galaxies in clusters (LP) is discrepant at more than the 2 — a level. 



to clusters based on other distance indicators will soon tell whether this early result is a 

^ 2(7 statistical fluke (Watkins & Feldman 1995; Strauss et al. 1995), whether the errors 
were underestimated, or whether something is systematically difi^erent between the BCG 
distances and the distances measured by other indicators. 

Shown in comparison are the expected rms bulk velocity in a standard CDM model 
(Q = 1, n = 1) normalized to COBE, for h — 0.5 (bottom) and 0.8 (top) (Sugiyama 1995). 
These theoretical curves would not change much if a 20% hot component is mixed with the 
cold dark matter, or if Vl is lower but still in the range 0.2 — 1.0, as long as n ^ 1. The main 
effect of Vt and Hq on P{k) is via kpeak- The predicted bulk velocity over ~ 100 /i^^Mpc 
is effectively an integral of P{k) over k < kpeak, and is therefore relatively insensitive to fl 
while it is quite sensitive to n. When compared to a theoretical prediction, the error should 
also include cosmic scatter due to the fact that only one sphere has been sampled from a 
random field. These errors are typically on the order of the measurement errors. 

The measurements of CMB fluctuations on scales < 90° are independent of the local 
streaming motions, but GI predicts an intimate relation between their amplitudes. The CMB 
fluctuations are associated with fluctuations in gravitational potential, velocity and density 
in the surface of last scattering at ^ ~ 10^, while similar fluctuations in our neighborhood 
have grown by gravity to produce the dynamical structure observed locally. The comparison 
between the two is therefore a crucial test for GL 

Before COBE, the local streaming velocities served to predict the expected level of 
CMB fluctuations. The local surveyed region of ~ 100 /i~^Mpc corresponds to a ~ 1° patch 
on the last-scattering surface. An important effect on scales > 1° is the Sachs- Wolfe effect 
(1967), where potential fluctuations induce temperature fluctuations via gravitational 
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redshift, ST/T = A^g/{3c'^). Since the velocity potential is proportional to $g in the linear 
and mildly-nonlinear regimes, A$g ~ Vx, where x is the scale over which the bulk velocity 
is V. Thus 5T/T > V^x/(3c^). A typical bulk velocity of ~300kms~^ across ~100/i~^Mpc 
(§ |3.1| ) corresponds to 6T/T> 10~^ at ~ 1°. If the fluctuations are roughly scale-invariant 
{n = 1), then 6T/T > 10^^ is expected on all scales > 1°. Bertschinger, Gorski & Dekel 
(1990) produced a crude 6T/T map of the local region as seen by a hypothetical distant 
observer, and predicted 6T/T> 10~^ from the local potential well associated with the GA. 
An uptodate version of the ST/T maps is provided by Zaroubi et al. (1997b), who added a 
proper treatment of the acoustic effects on sub-degree scales for various cosmological models. 

Now that CMB fluctuations of ~ 10~^ have been detected practically on all the relevant 
angular scales, the argument can be reversed: if one assumes GI, then the expected bulk 
velocity in the surveyed volume is ~300kms~^, i.e., the inferred motions of § @ are most 
likely real. If, alternatively, one accepts the peculiar velocities as real for other reasons, then 
their consistency with the CMB fluctuations is a relatively sensitive and robust test of the 
validity of GI. This test is unique in the sense that it addresses the specific fluctuation growth 
rate as predicted by GI theory (§ ISTTp . It is robust in the sense that it is quite insensitive to 
the values of the cosmological parameters and is independent of the complex issues involved 
in the process of galaxy formation. 



3.2. Power Spectrum from the Velocity Field via POTENT 

One way to compute the power spectrum is via the smoothed mass density field as 
recovered by POTENT (Kolatt & Dekel 1997). The key is to correct the result for systematic 
deviations from the true P{k). The data suffers from distance errors and sparse, nonuniform 
sampling, and they were heavily smoothed. The P{k) is computed from within a window of 
effective radius ~ 50 h~^Mpc, say, where the densities are weighted inversely by the squares 
of the local errors. The density field is zero-padded in a larger periodic box in order to 
enable an FFT procedure. The P{k) is computed by averaging the amplitudes of the Fourier 
transforms in bins of k. This procedure yields an "observed" P{k), which we term 0{k). 

The systematic errors in the above procedure are then modeled by 0{k) = M{k)[S{k) + 
N{k)], where S{k) is the true signal P{k), N{k) is the noise, and M{k) represents the effects 
of sampling, smoothing, applying a window etc. The correction functions M{k) and N{k) 
can be derived from Monte Carlo mock catalogs (§ |2.6| ). The factor M{k) is derived first by 
= S/ {O)no-noise, whcrc S here is the known power spectrum built into the simulations, 
and the averaging is over mock catalogs not perturbed by noise. Then N{k) is computed by 
N = M~^{0) noise — S, where the averaging is over noisy mock catalogs. 

Equipped with the correction functions M{k) and N{k), the P{k) observed from the 
real universe, 0{k), is corrected to yield the true P{k) by S{k) = M{k)^^0{k) — N{k). The 
recovered mass-density P{k) is shown in Figure ^ in three thick logarithmic bins covering the 
range 0.04 < k < 0.2 ( /i~^Mpc)~^, within which the results are reliable. The robust result is 
P{k)f{nf = (4.6±1.4) X 10^( /i^^Mpc)^ atk = 0.1( /i-^Mpc)-^ (using the convention where 
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Fig. 6. — The mass power spectrum [x/(r2)^] from POTENT Mark III velocities (filled symbols), 
with la random errors. The curves are COBE-normalized theoretical predictions for flat CDM 
with h = 0.5 (solid) and h = 0.7 (dashed). Left: A-CDM with n = 1 and Q growing upwards. 
Right: tilted-CDM with = 1 and n growing upwards. Typical cosmic scatter (CS) for = 1 and 
n = 1 is indicated (Kolatt & Dekel 1997). 



the Fourier transform is defined with no 27r factors in it's coefficient). The logarithmic slope 
at k = 0.1 is — 1.45±0.5. This translates to as ^^'^ — 0.7 — 0.8, depending on where the peak 
in P{k) is (see § p]). 

The observed P{k) is compared in the figure to the linear predictions of a family of 
Inflation-motivated flat CDM models {Q + = 1), normalized by the 4-year COBE data, 
with the Hubble constant fixed at h = 0.5 or h = 0.7. For n = 1, maximum likelihood is 
obtained at ~ O.Th^Q-^ ±0.1. For Q = 1, assuming no tensor fluctuations, the linear power 
index is n ~ 0.75h7^-^ ± 0.1. 



3.3. Power Spectrum from Velocities for COBE-CDM Models 

The power spectrum, in a parametric form including Q, h and n among the parameters, 
has alternatively been determined from the velocity data via a Baysian likelihood analysis 
(Zaroubi et al. 1997a; see also Kaiser & Jaffe 1995). According to Bayse, the probability of 
the model parameters (m) given the data (d), which is the function one wishes to maximize, 
can be expressed as P{m\d) = P{d\rn) P{m) / P{d). The probability P{d) serves here as a 
normalization constant. Without any external constraints on the model parameters, one 
assumes that P{m) is a constant in a given range. The remaining task is to maximize the 
likelihood C = P{d\m) as a function of the model parameters. This function can be written 
down explicitly. 

Under the assumption that the velocities and the errors are both Gaussian random fields 
with no mutual correlations, the likelihood can be written as £ = {2n\D\)-'^/'^ exp{-diD-j^dj/2), 
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where di are the data at points i = 1, ...,N, and Dij is the covariance matrix, which can be 
spht into covariance of signal (s) and covariance of noise (n), Dij = {didj) = {siSj) + {riirij). 
If the errors are uncorrelated, the noise matrix is diagonal. The signal matrix is computed 
from the model P{k)a.s a function of the model parameters. 

Zaroubi et al. (1997a) used a parametric model for the PS of the general form Pk = 
A k"- T [r i] k) , where T{k) is a small-scale filter of an assumed shape characterized by free 
parameters Fj, fc" is the initial P{k) which is still valid on large scales today, and A is a 
normalization factor. The normalization can either be determined by COBE's data (for 
given Q, A, h, n and tensor/scalar fluctuations), or be left as a free parameter to be fixed by 
the velocity data alone. The filter T{k) can either be taken from a specific physical model 
{e.g., CDM, where F = Qh), or be an arbitrary function with enough flexibility to fit the 
data. 

The robust result for all the models is a relatively high amplitude, with P{k)f{Q)'^ = 
(4.8 ± 1.5) X 10^( h~^Mpc)^ at = 0.1( h~^Mpc)~^ . An extrapolation to smaller scales using 
the different CDM models yields ag^^'^ = 0.88 ± 0.15 (for the dispersion in top-hat spheres 
of radius 8 h~^Mpc). 

Within the general family of CDM models, allowing for a cosmological constant in a fiat 
universe and a tilt in the spectrum, the parameters are confined by a 90% likelihood contour 
of the sort Q h^^^ n'^ = 0.8 ±0.2, where /i = 1.3 and u = 3.4, 2.0 for models with and without 
tensor fluctuations respectively. Figure |^ displays the likelihood map in the Q — n plane for 
these models. For open CDM the powers are /i = 0.95 and u = 1.4 (no tensor fluctuations). 
A F-shape model free of COBE normalization yields only a weak constraint: F = 0.4 ± 0.2 
(where F is not necessarily Qh). 




Fig. 7. — Likelihood contour maps in the — n plane for flat CDM with and without tensor 
fluctuations and h = 0.65. Contour spacing is unity in log- likelihood. Under the assumption of 
a distribution, the two-dimensional 90 percentile corresponds to 2.3 contours, and the one- 
dimensional 90 percentile corresponds to 1.35 contours (Zaroubi et al. 1997a). 



Both Q and n obtained by the likelihood analysis from the raw peculiar velocities tend 
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to be slightly higher (~ 20%) than their estimates based on the P{k) recovered from the 
POTENT output. This difference may arise from the different relative weighting assigned to 
the different wavelengths in the two analyses. The difference between the results obtained in 
the two different ways is on the order of the errors in each analysis and the cosmic scatter. 
Very similar estimates of P{k) are obtained from a preliminary analysis of the SFI sample 
of Sc galaxies (in preparation). 

In summary: The "standard" CDM model is marginally rejected at the ~2(j level, while 
each of the following modifications lead to a good fit to the peculiar velocities and large-scale 
CMB data: n ^ 1, Qi, ^ 0.3, or f2 ^ 1. The strong implication on the dark matter issue 
is that values of fl as low as ~ 0.2 are ruled out with high confidence (independent of A), 
leaving, in particular, no room for the baryonic PIB model. 



3.4. Peculiar Velocities vs Small-scale CMB Fluctuations 

Sub-degree angular scales at the last scattering surface correspond to the < 100 h~^Mpc 
comoving scales explored by peculiar velocities today. Thus, under the assumption that 
the local neighborhood is typical, the power spectrum on these scales is simultaneously 
constrained by the mass-density fluctuations in our cosmological neighborhood and by the 
CMB fluctuations. 

The sub-degree CMB fluctuations are now being explored by many balloon-born exper- 
iments, and in less than a decade we expect accurate results from the CMB satellites MAP 
and Planck. These measurements will eventually allow a simultaneous likelihood analysis 
of the two kinds of data. At this point, however, although there are already preliminary 
detections of the first acoustic peak in the angular power spectrum, the uncertainties are 
still large. Any current comparison is therefore limited to the semi-quantitative level. The 
range of parameters permitted by the peculiar velocity data for power spectra of the CDM 
family (Zaroubi et al. 1997a) can be translated to a range of angular power spectra, Q. This 
range is plotted against current observations in Figure |^. 

The immediate conclusions from a visual inspection of the figure are that a wide range 
of CDM models can simultaneously obey the two data sets, but there is a subset of models 
that fits the velocities well but seems to fail to produce a high enough acoustic peak in the 
CMB spectrum. 

The acoustic peak in the CMB is sensitive to Vt^ and the observations prefer a high 
baryon content, ~ 0.025 (similar to the value measured by Tytler et al. 1996; Buries & 
Tytler 1996), while the peculiar velocities have little to add because P{k) is hardly affected 
by Vt^. 

The power index n is important in both cases; the peculiar velocities allow values of n 
significantly lower than unity, but the current CMB data seem not to tolerate values of n 
below 0.9 or so. 

The peculiar velocity data prefers VL > 0.4, and the location of the first acoustic peak 
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Fig. 8. — Angular power spectrum in the CMB. The current data (symbols) are compared to 
certain CDM models that fit the peculiar velocity data at 90% likelihood. The extreme models 
shown are either with low h and high Q,i,h? or vice versa. Left: tilted CDM with = 1, including 
tensor fluctuations. Right: Flat CDM with n = 1 (Zaroubi et al. 1997b). xxxxxx 



in the sub-degree CMB data indicates in agreement a high value of i7 + ^l^- 



4. DIRECT MEASUREMENTS OF Q FROM PECULIAR VELOCITIES 

Assuming that the inferred motions are real and generated by GI, they can be used to 
estimate Q in several different ways. Most of the evidence from virialized systems on scales 
< 10 h^^Mpc suggest a low mean density of fi~0.2 (see Dekel, Burstein & White 1997). The 
spatial variations of the large-scale velocity field provide ways to measure the mass density 
in a larger volume that may be closer to a "fair" sample. One family of such methods is 
based on comparing the dynamical fields derived from velocities to the fields derived from 
galaxy redshifts (§ ^). These methods can be applied in the linear regime but they always 
rely on the assumed biasing relation between galaxies and mass often parameterized by b, 
so they actually provide an estimate of f3 = f{Q)/b. Another family of methods measures (3 
from redshift surveys alone, based on 2r-space deviations from isotropy (see Strauss 1997). In 
the present section, we focus first on methods that rely on non-linear effects in the peculiar 
velocity data alone, and they thus provide estimates of Q independent of galaxy density 
biasing. These methods are based on the assumption that the initial fluctuations were 
Gaussian. 



4.1. Divergence in Voids 

A diverging flow in an extended low-density region can provide a robust dynamical lower 
bound on fl, based on the fact that large gravitating outflows are not expected in a low-Jl 
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universe (Dekel & Rees 1994). In practice, for any assumed value of Q, the partial derivatives 
of the smoothed observed velocity field are used to infer a non-linear approximation for the 
mass density via the approximation • A key point is that this approximation is typically 
an overestimate, 6c > S (when the true value of fl is used). For fluctuations that started 
Gaussian, the probability that 6c is an overestimate, in the range 6 < —0.5, is well over 99%. 
Analogously to (5o ~ — ^^^'^V-w, the 6c inferred from a given diverging velocity field is more 
negative when a smaller Q is assumed, so for a small enough Q one may obtain 5 < — 1 in 
certain void regions. Such values of 6 are forbidden because mass is never negative, so this 
provides a lower bound on Q. 

The inferred 6c field, smoothed at 12/i^^Mpc, and the associated error field as, were 
derived by POTENT at grid points from the observed radial velocities of Mark III. Focusing 
on the deepest density wells, the input Q was lowered until 6c became significantly smaller 
than —1. The most promising "test case" provided by the Mark III data is a broad diverging 
region centered near the supergalactic plane at the vicinity of (X, Y) = (—25, —40) in h~^Mpc 
— the "Sculptor void" of galaxies between the GA and the "Southern Wall" extension of 
PP (Figure |^, compare to Fig. ^). 




X (1000 km/sec) X (1000 km/sec) 

Fig. 9. — Maps of 6c inferred from the observed velocities near the Sculptor void in the Super- 
galactic plane, for two values of 0. The LG is marked by and the void is confined by the Pavo 
part of the GA (left) and the Aquarius extension of PP (right). Contour spacing is 0.5, with 6c = 
heavy, 6c > solid, and 6c < dotted. The heavy-dashed contours mark the illegitimate downward 
deviation of 5c below —1 in units of as, starting from zero (i.e., 5c = —1) and decreasing with 
spacing — O.So"^. The value O = 0.2 is ruled out at the 2.9a level (Dekel & Rees 1994). 

Values of r2^1 are perfectly consistent with the data, but 6c becomes smaller than — 1 
already for Q = 0.6. The values fi = 0.3 and 0.2 are ruled out at the 2.4cr and 2.9a levels in 
terms of the random error as. 

This result is still to be improved. The systematic errors have been partially corrected 
in POTENT, but a more specific investigation of the biases affecting the smoothed velocity 
field in density wells is still in progress. For the method to be effective one needs to find a 
void that is (a) bigger than the correlation length for its vicinity to represent the universal 
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Q, (b) deep enough for the lower bound to be tight, (c) nearby enough for the distance errors 
to be small, and (d) properly sampled to trace the velocity field in its vicinity. 

Note that this method does not require that the void be spherical or of any other 
particular shape, and is independent of galaxy density biasing. Another pro is that there is 
no much cosmic scatter — one deep and properly sampled void is enough for a meaningful 
constraint. The main limitation is the poor (and perhaps biased) sampling of the velocity 
field in the vicinity of a void. 



4.2. Deviations from Gaussian PDF 

Assuming that the initial fluctuations are a Gaussian random field, the one-point prob- 
ability distribution function (PDF) of smoothed density develops a characteristic skewness 
due to non-linear effects early in the non-linear regime {e.g., Kofman et al. 1994). The 
skewness of 5 is given according to second-order perturbation theory by 

(53)/(52)2^(34/7-3-n) , (9) 

with n the effective power index of the power spectrum near the (top-hat) smoothing scale 
(Bouchet et al. 1992). Since this ratio of moments for 5 is practically independent of f2, and 
since V-v~— the corresponding ratio for V-v must strongly depend on f2, and indeed in 
second-order it is (Bernardeau et al. 1995) 

T3 ^ {{V-vf)/{{V-v)Y ^ -fin)-\26/7 - 3 - n) . (10) 

Using N-body simulations and 12/i~^Mpc smoothing one indeed finds T3 = —1.8 ± 0.7 for 
Q=l and T3 = —4.1 ±1.3 for = 0.3, where the error is the cosmic scatter for a sphere of 
radius AOh'^Mpc in a CDM universe {Hq = 75, 6=1). An estimate of T3 in the current 
POTENT velocity field within 40 h~^Mpc is —1.1 ±0.8, where the errors this time represent 
distance errors. With the two errors added in quadrature, fl = 0.3 is rejected at the ~ 2a 
level (somewhat sensitive to the assumed P{k)). 

Since the present PDF contains only part of the information stored in the data and is 
in some cases not that sensitive to the initial PDF (IPDF), a more powerful bound can be 
obtained by using the detailed present velocity field v{x) to recover the IPDF, and using 
the latter to constrain Q by measuring it's f2-dependent deviation from the assumed normal 
distribution (Nusser & Dekel 1993). The necessary "time machine" is provided by the 
Eulerian interpretation of the Zel'dovich approximation (Nusser & Dekel 1992). 

The velocity out of POTENT Mark II, within a conservatively selected volume, was 
fed into the IPDF recovery procedure with Q either 1 or 0.3, and the errors due to distance 



errors and cosmic scatter were estimated. Figure |T0| shows the recovered IPDF's. The IPDF 
recovered for = 1 is marginally consistent with Gaussian, while the one recovered for 
Q = 0.3 shows significant deviations. The largest deviation, bin by bin in the IPDF, is ^2a 
for Q = l and >4a for Q = 0.3, and a similar rejection of f2 = 0.3 is obtained with a x^-type 
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Fig. 10.— Q from IPDF. The density IPDF recovered from tlie G2 POTENT peculiar velocity 
field (solid) , compared to a normal distribution (short dash) and with the IPDF recovered from the 
velocity field of Gaussian CDM simulations (triangles). The assumed VL is 1.0 (left) or 0.3 (right). 
The simulations are of fio = accordingly (Nusser k. Dekel 1993). 



statistic. The skewness and kurtosis are poorly determined because of noisy tails but the 
replacements and allow a rejection of = 0.3 at the (5 — 6)0" levels. 

The main advantage of the methods based on the PDF is their insensitivity to galaxy 
density biasing. The main weakness is the need for a "fair" sample; the cosmic scatter is 
large due to the large smoothing scale within the limited volume. 



5. GALAXY DENSITY VS. VELOCITIES: O AND BIASING 
5.1. Galaxies vs Mass: Fit of GI and Linear Biasing Model 

The theory of GI combined with the assumption of linear biasing for galaxies predict a 
correlation between the dynamical density field and the galaxy density field, which can be 
addressed quantitatively based on the mock catalogs and the estimated errors in the two data 
sets. Figure |Tl] compares density maps in the Supergalactic plane for IRAS 1.2 Jy galaxies 
(5g) and POTENT Mark III mass (5), both G12 smoothed. The general correlation is evident 
— the GA, PP, Coma and the voids all exist both as dynamical entities and as structures of 
galaxies. To evaluate goodness of fit, Figure P shows the statistic = N'^ {Sc-bS)^ /a^ 



as computed from the data in comparison with its distribution over pairs of Mark III and 
IRAS 1.2 Jy mock catalogs. The fact that the data lies near the center of this distribution 
indicates that the two data sets are consistent with being noisy versions of an underlying 
fluctuation field and that the data are in agreement with the hypotheses of GI plus linear 



biasing (Dekel et al. 1993; Sigad et al. 1997; more in § |5^) . 



What is it exactly that one can learn from the observed v — 5q correlation (Babul et 
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Fig. 11. — Mass versus galaxies. POTENT mass (17 = 1) versus IRAS galaxy density fields in the 
Supergalactic plane, both smoothed G12. Contour spacing is 0.2. The heavy contour marks the 
boundary of the comparison volume of effective radius 46/i~^Mpc. The height in the surface plot 
is proportional to 5. The LG is at the center, GA on the left, PP on the right, and the Sculptor 
void in between (Sigad et al. 1997). 

al. 1994)? First, it argues that the velocities are real because it is hard to invoke any other 
reasonable way to make the galaxy distribution and the TF measurements agree so well. 
On the other hand, although it is true that gravity is the only long-range force that could 
attract galaxies to stream toward density concentrations, the fact that eiV — 5o correlation is 
predicted by GI plus linear biasing does not necessarily mean that the observation can serve 
as a sensitive test for either. Recall that converging (or diverging) flows tend to generate 
overdensities (or underdensities) simply as a result of mass conservation, independent of the 
source of the motions. 

Let us assume for a moment that galaxies trace mass, i.e., that the linearized continuity 
equation, 5 — —V-V, is valid for the galaxies as well. The observed correlation (in the hnear 
approximation) is then S oc — V-u, and together they imply that SocS, or equivalently that V-t; 
is proportional to its time average. This property is not exclusive to GI; one can construct 
counterexamples where the velocities are produced by a non-GI impulse. 
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Fig. 12. — Goodness of fit in the comparison of G12 density fields of POTENT mass and IRAS 
galaxies, as tested by a statistic that is an error weighted sum of differences between the two 
fields (Sigad et al. 1997). 



Even irrotationality does not follow from 5oc— V-v; it has to be adopted based on the- 
oretical arguments in order to enable reconstruction from radial velocities or from observed 
densities. Once continuity and irrotationality are assumed, the observed 5oc— V-t> implies 
a system of equations which is identical in all its spatial properties to the equations of GI, 
but can differ in the constants of proportionality and their temporal behavior. It is there- 
fore impossible to distinguish between GI and a non-GI model which obeys continuity plus 
irrotationality based only on snapshots of present-day linear fluctuation fields. This makes 
the relation between CMB fluctuations and velocities an especially important test for GI. 

On the other hand, the fact that the constant of proportionality in 5 oc — V-u is indeed the 
same everywhere is a non-trivial requirement from a non-GI model. For example, a version 
of the explosion scenario (Ostriker & Cowie 1981; Ikeuchi 1981), which tested successfully 
both for irrotationality and v — 5 correlation, requires special synchronization among the 
explosions (Babul et al. 1994). 

So, what is the v—5g relation good for? While it's sensitivity to GI is only partial, this 
relation turns out to be quite sensitive to the validity of a continuity-like relation for the 
galaxies. When the latter is strongly violated all bets are off for the v — Sq relation. A non- 
linear biasing scheme would make continuity invalid for the galaxies, which would ruin the 
v—Sq relation even if GI is valid. The observed correlation is thus a sensitive test for density 
biasing. It implies, subject to the errors, that the ~ 12 /i~^Mpc-smoothed density fields of 
galaxies and mass are related via a biasing relation that could be crudely approximated by 
a linear relation with b of order unity (but see a refinement of this in § |5.5| ). 

Now that the data of peculiar velocities and the data of galaxy density are found to be 
compatible with the model of GI and linear biasing, they can be combined to constrain the 
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degenerate parameter j3 = Vf''^ /h. The comparison between the two data sets can be done 
in several different ways. In particular, it could be done by comparing density fields derived 
locally from the two data sets (§ [5.2| ), or by comparing velocities derived from the two data 
sets (§ p.4| ). It can be done successively by first recovering fields from each data set and then 
combining them to obtain (3, or by a simultaneous recovery of fields and beta determination 
from the two data sets (§ |5.3| ). It can be done by direct comparison of fields in r-space, or 
by comparing coefficients mode by mode in a model expansion in z-space (Davis, Nusser & 
Wilhck 1996). 



5.2. Density-Density Comparison on Large Scales: POTIRAS 

The main advantage of comparing densities is that they are local. The densities are 
independent of long-range effects due to the unknown mass distribution outside the sampled 
volume, which could affect the velocities. The densities are also independent of reference 
frame, and can be reasonably corrected for non-linear effects. 

The POTENT analysis extracts from the peculiar velocity data a mildly-nonlinear mass 
density fluctuation fleld in a spatial grid, smoothed G12 (§ [ZS] ). The associated real-space 
density fleld of galaxies can be extracted with similar smoothing from a whole-sky redshift 
survey such as the IRAS 1.2 Jy survey (see Strauss & Willick 1995; Sigad et al. 1997). 

A brief summary of the recovery of the IRAS density fleld is as follows. The solution to 
the linearized GI equation V-v = — /5 for an irrotational fleld is 

v{x) = ^ I d^x'6{x') ,^''^,^ . (11) 

47r Jail space \x' - X]'^ 

The velocity is proportional to the gravitational acceleration, which ideally requires full 
knowledge of the distribution of mass in space. In practice, one is provided with a flux- 
limited, discrete redshift survey, obeying some radial selection function 0(r). The galaxy 
density is estimated by l+6G{x) = J2'n~^(f){ri)~^6%^a^{x — Xi), where n = V^^ l]0(n)~^ is the 
mean galaxy density, and the inverse weighting by restores the equal-volume weighting. 
Eq. |ll| is then replaced by 

Under the assumption of linear biasing, the cosmological dependence enters through (3. The 
integration is limited to r < Rmax where the signal dominates over shot-noise. S{y) is a 
small-scale smoothing window (>500kms~^) essential for reducing the effects of non-linear 
gravity, shot-noise, distance uncertainty, and triple-value zones. 

The distances are estimated from the redshifts in the LG frame by 



/•, = Zi - it, ■ \v{x^ - i)(0)] 



(13) 
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Equations 0-|T3| can be solved iteratively: make a first guess for the Xj, compute the Vi 
by Eq. 0, correct the Xi by Eq. |T3|, and so on until convergence. The convergence can be 
improved by increasing f3 gradually during the iterations from zero to its desired value. 

Even under 12 h~^Mpc smoothing, Sq is of order unity in places, necessitating a mildly- 
nonlinear treatment. Local approximations from v to 6 were discussed in § |2.5| , but the 
non-local nature of the inverse problem makes it less straightforward. A possible solution 
is to find an inverse relation of the sort V-t? = F{Q,6g), including non-linear gravity and 
non-linear biasing. This is a Poisson-like equation in which —i36g{x) is replaced by F{x), 
and since the smoothed velocity field is still irrotational for mildly-nonlinear perturbations, 
it can be integrated analogously to Eq. |ll|. With smoothing of 10 /i~^Mpc and (3 = 1, the 
approximation based on an empirical inverse to 6c has an rms error < 50kms^^. 

In recent applications, the galaxy density field is recovered from the noisy IRAS data via 
a Power- Preserving Filter (PPF, by A. Yahil, described in Sigad et al. 1997) — a modification 
of the Wiener Filter. The PPF returns a field that is not far from the Wiener, most probable 
field, but it makes the result more realistic by forcing the variance to be constant in space 
despite the fact that the errors vary. 




Fig. 13. — /3 from POTENT vs IRAS density comparison. The smoothing is G12 and the compar- 
ison volume is of effective radius 40/i~^Mpc. Two-dimensional regression lines are marked. Left: 
Regression between the averages of 20 mock catalogs of each type, showing a bias as small as 4%. 
Right: Real data (Sigad et al. 1997). 



The simplest way of comparing the POTENT and IRAS density fields is via a two- 
dimensional linear regression using the values of the fields at grid points within a local 
comparison volume. The errors of both fields enter the regression. The comparison volume 
is determined by equal-error contours. 

The latest comparison of POTENT Mark III and IRAS 1.2 Jy data at 12 /i^^Mpc 
smoothing within a volume of (65 /i~^Mpc)'^ yields /57 = 0.86 ± 0.12 (Sigad et al. 1997). The 



corresponding scatter diagram is shown in Figure 13. The systematic error in this derivation 
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of only 4%, is deduced from the analogous scatter diagram for the averages of 20 random 
mock catalogs of each type. This is an update of the higher estimate Pj = 1.3 ± 0.3 (Dekel et 
al. 1993) obtained based on an earlier version of POTENT with the Mark II velocities and 
the IRAS 1.9 Jy redshifts. 

Similar comparisons of the mass density field with the density of optical galaxies indicate 
a similar correlation and a ~30% lower estimate for (3o (Hudson et al. 1995), in agreement 
with the ratio of biasing factors, ho/h^ 1.3, obtained by direct comparison of optical and 
IRAS galaxy densities. 

A comparison of similar nature of the POTENT Mark III data with the density dis- 
tribution of Abell/ACO i? > clusters and the corresponding predicted velocities at G15 
smoothing yields similar consistency out to distances of ~ 60 /i~^Mpc, and an estimate of 
(3c = 0.26 ± 0.11 (Plionis et al. 1997). This is consistent with a linear biasing factor for the 
clusters that is about 4 times larger than that of galaxies, in accordance with the observed 
ratio of about 4^ for the corresponding correlation functions (see Bahcall 1997). 
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Fig. 14. — /3 from power spectra of galaxies versus mass. The estimates from various galaxy 
density samples were all translated from redshift to real space using Kaiser's approximation and 
the best-fit value of (5. The P{k) from peculiar velocities (Fig. ^) is marked by solid symbols for 
= 1 (open symbols for = 0.3). The values of fi (for any Q) can be read directly from the 
vertical offset of the solid symbols and the corresponding curves (Kolatt & Dekel 1997). 



A direct comparison of the mass power spectrum as derived from peculiar velocities 
]2]) with the galaxy power spectra as derived from different redshift and angular surveys 
is shown in Figure [l^ (Kolatt & Dekel 1997). It demonstrates a similarity in shape and 
yields for the various galaxy types (3 values in the range 0.77 — 1.21, with a typical error of 
±0.1. For IRAS galaxies typically /3 ^ 1, and for optical galaxies /? ^ 1. These estimates 
do not directly address the value of f2, but it is clear from the figure that if Vt is as small as 
~0.3, then all galaxy types must be severely antibiased. 

In principle, the degeneracy of Vt and h is broken in the mildly-nonlinear regime, where 
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S{v) is no longer oc /~^. Compatible mildly-nonlinear corrections in POTENT and in the 
IRAS analysis allowed a preliminary attempt to separate these parameters using the Mark II 
and IRAS 1.9 Jy data (Dekel et al. 1993). Unfortunately, nonlinear biasing effects are hard 
to distinguish from non-linear gravitational effects, so a specific nonlinear biasing scheme is 
a prerequisite for such an analysis. 



5.3. Simultaneous Fit of Velocity and Density: SIMPOT 

The dynamical fields and (3 can be recovered simultaneously by a fit of a parametric 
model for the potential field to the combined data of the observed radial peculiar velocities 
and the distribution of galaxies in redshift space. This procedure takes advantage of the 
complementary features of the data in the recovery of the fields, it enforces the same effective 
smoothing on the data without preliminary reconstruction procedures such as POTENT, and 
it obtains a more reliable best fit by simultaneous rather than successive minimization. It 
has been implemented so far to the forward TF data, but it can be generalized in principle 
to minimize inverse TF residuals. 

In the SIMPOT procedure by Nusser & Dekel (1997), the model for the potential field 
is taken to be an expansion in spherical harmonics Yi^ and Bessel functions ji where the 
coefficients ^Imn are the free parameters. 



*W = E E E ^Irnnjl 



(14) 



;=1 m=—l n=l 



The model radial velocity is derived from this potential hj u = —d^/dr, and the model 
density in redshift space is derived using linear theory via 5 = f^^V^^+s~'^{d/ds){s'^d^/ds). 
The second term refiects redshift distortions, where s is the radial variable in redshift space. 
The resulting models for u{x) and 6{s) are expansions in certain functions Aimn and Bij^n 
that arc appropriate combinations of the original base functions. 

The combined to minimize as a function of the parameters ^imn and /3 is the sum of 



and 



/ 
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(15) 



(16) 



The observations arc the peculiar velocities uf'^ and a continuous density field in redshift 
space 5°^'^{s)^ that is somewhat more tricky to obtain. The j3 dependence enters only in x^, 
both via the velocity-density relation and the redshift distortions. 

A SIMPOT fit to the Mark III peculiar velocity data and the IRAS 1.2 Jy redshift survey 
provides first hints for scale-dependent biasing: /3jf=i0.6 and 1.0 (±0.1) for smoothings that 
are roughly equivalent to Gaussian with radii 6 and 12 /i~^Mpc respectively. 
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5.4. Velocity- Velocity Comparison on Small Scales: VELMOD 

Earlier comparisons of the peculiar velocities from the Mark II catalog and the velocities 
predicted from the IRAS redshift surveys (QDOT, 1.9 Jy, and 1.2 Jy) yielded estimates of 
/3i in the range 0.4 - 1.0 (Kaiser et al 1991; Nusser & Davis 1994). 

The more sophisticated recent VELMOD method of comparison (Willick et al. 1997b) 
compares the raw peculiar velocity data with a "model" velocity field that is predicted from 
the IRAS 1.2 Jy redshift survey. It is done without attempting to reconstruct a velocity field 
from the data. The key feature of VELMOD is that it explicitly allows for a non-unique 
mapping between real space and redshift space. Triple valuedness in the redshift field as 
well as non-negligible small-scale velocity "temperature" are treated in a unified way. This 
is done by expressing the probability that an object at distance r has a redshift z by 



P{z\r) = — exp 

V 27ro",, 



1 [z — r — u{r)]' 
'2 2^! 



(17) 



where u{r) is the radial component of the model velocity field and is the small-scale 
velocity noise (which can in principle be a function of position). The above probability is 
then multiplied by the TF probability factor, P(m, t], r), and integrated over the entire line- 
of-sight to obtain the probability of the observable quantities {m,7], z). One then maximizes 
that probability over the entire data set. 

The method is computer-intensive because numerical integrals are required for each 
galaxy, and for each fit parameter (TF parameters, au, velocity field model parameters, 
etc.). This effort is worthwhile to the degree that the velocity field is triple- valued or the 
small-scale noise is comparable to the TF error. In particular, VELMOD is more rigorous 
in an analysis of the very local {z < 3000 km s~^) region. 

VELMOD has been applied with a Gaussian smoothing of 3/i~^Mpc to the IRAS 1.2 
Jy redshift survey and a subset of 838 spiral galaxies from the Mark III catalog within 
z < 3000 km s~^ of the Local Group. The method was tested successfully using mock cat- 
alogs drawn from the N-body simulation of Kolatt et al. (1996). When applied to the real 
data it yielded consistency with the model of linear GI and linear biasing once an artificial 
quadrupole was allowed, with (5j = 0.5 ± 0.1 at 3 /i~^Mpc. The catch is that it is not at 
all clear why linear GI and the simplified deterministic biasing should be valid for the den- 
sities and velocities at such high resolution. The estimated value of (3 should therefore be 
interpreted with caution. 

It is interesting to note that Shaya et al. (1995) obtained a similarly low value for (3o 
from the same local neighborhood. They applied the least-action reconstruction method to 
a redshift survey of several hundred spirals in comparison with TF data. Their method is 
likely to underestimate (3 because it assumes that the mass is all concentrated in the centers 
of galaxies and groups and thus tends to overestimate the gravitational forces between them. 
This systematic effect is yet to be quantified. 



Table |^ summarizes the estimates of /3 and fl from cosmic flows. 
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Table 1: O and /3 from Cosmic Flows 



Peculiar Gaussian IPDF Nusser & Dekel 93 Q > 0.3 f> 4(t) 

Velocities Skcwness(V-u) Bernardeau et al. 94 Q > 0.3 (2(t) 

Alone Void Dekel & Rees 94 n> 0.3 (2.4ct) 

Power spectrum Kolatt & Dekel 97 as^^'^ = 0.7 ± 0.15 

+COBE Zaroubi et al. 97a as^^'^ = 0.8 ± 0.15 



Galaxy 
Density 
vs. 

Velocities 



M2-QD0T V 
M3-I1.2 v-dipole 
M3-I1.2 v-inverse 
M3-I1.2 V G3 
M3-I1.2 S G12 
M3-I1.2 5/v G6-12 

M2-0ptical V 
TF-Optical 
M3-0ptical S G12 



Kaiser et al. 91 
Nusser & Davis 94 
Davis et al. 96 
Willick et al. 96 
Sigad et al. 97 
Nusser & Dekel 96 

Hudson 94 
Shaya et al. 94 
Hudson et al. 95 



Pi = o.9l[l:?5 

= 0.6 ± 0.2 
Pj = 0.6±0.2(?) 
pj = 0.5 ±0.1 
pj = 0.86 ±0.15 
/3, = 0.6 - 1.0 scale 



Po 
Po 
Po 



0.5 ±0.1 
0.35 ±0.1 
0.75 ± 0.2 



M3-clusters G15 



Plionis et al. 97 



Pc 



0.26 ±0.11 



Redshift 


C 11.2 


Peacock & Dodds 94 




= 1.0 ±0.2 


Distortions 


C 11.2 


Fisher et al. 94a 


Pi 


= 0.45l°i 




Ylm 11-2 


Fisher et al. 94b 


Pi 


= 1.0 ±0.3 




Pk 11.2. QDOT 


Cole et al. 95 


Pi 


= 0.5 ±0.15 




{ 11.2, QDOT 


Hamilton 95 


Pi 


= 0.7 ±0.2 




Ylm 11.2 


Heavens & Taylor 95 


Pi 


= 1.1 ±0.3 




Pk 11.2 


Fisher & Nusser 96 


Pi 


= 0.6 ±0.2 



CMB 
Dipole 


vs galaxies angular 
vs galaxies redshift 


Yahil et al. 86 
Strauss et al. 92 
Rowan- Rob. et al. 91 


Pi 
Pi 
Pi 


= 0.9 ±0.2 
= 0.4-0.85 

— n 8+0-2 

- U.»_o.l5 




vs galaxies angular 
vs galaxies redshift 


Lynden-Bell et al. 89 
Hudson 93 


Po 
Po 


= 0.3-0.5 

- 7+0-"^ 

— U. ( -0.2 




clusters 


Scaramella et al. 91 
Phonis et al. 91 


Pc 
Pc 


~ 0.13 

~ 0.17-0.22 



p = n^-yb, he : ho- hi ^ 4.5 : 1.3 : 1.0, a^VlP-^ = (0.69 ± 0.05)/?,, 

M3= Mark HI, 11.2 = IRAS 1.2 Jy, G12 = Gaussian smoothing 12 h~^Mpc 
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5.5. Galaxy Biasing as a Stochastic Process 



In all the methods described in § |], the cosmological parameter of interest Q is con- 
taminated by the uncertain relation between galaxy and mass density, the so called "galaxy 
biasing". Nontrivial galaxy biasing clearly exists. The fact that galaxies of different types 
cluster differently (Dressier 1980) implies that at least some do not trace the underlying 
mass. This is hardly surprising because any reasonable physical theory would predict non- 
trivial biasing (Kaiser 1984; Davis et al. 1985; Bardeen et al. 1986; Dekel & Silk 1986; Dekel 
& Rees 1987; Braun, Dekel & Shapiro 1988; Weinberg 1995). In particular, simulations of 
galaxy formation in a cosmological context {e.g., Cen & Ostriker 1992; 1993; Lemson et al. 
1997) indicate a biasing relation that is non-linear in density, is varying with scale, and has 
a statistical scatter reflecting dependencies on factors other than density. 

One should therefore not be surprised by the fact that the various estimates of /? span 
a large range, from less than one half to more than unity. Some of this scatter is due to the 
different types of galaxies involved, and some may be due to remaining effects of non-linear 
gravity or other systematic errors, but a significant fraction of the scatter in (3 is likely to 
reflect non-trivial properties of the biasing scheme. This means that translating a measured 
j3 into VL is non-trivial; it requires a detailed knowledge of the relevant biasing scheme. 

In order to strengthen this point, we demonstrate below that an obvious source of 
systematic variations in (3 is the inevitable statistical scatter in the biasing process (Dekel & 
Lahav 1997). This scatter in the relation between densities can be interpreted as reflecting 
the dependence of galaxy formation efficiency, or galaxy density, on physical properties 
of the protogalaxy environment other than density. These could be local properties such 
as the potential field, the deformation tensor, tidal effects, and angular momentum, or 
long-range effects carried by radiation or particles from neighboring sources. In the simple 
example below, we assume that this scatter in the biasing is local and neglect possible spatial 
correlations. 

Let Six) be the field of mass-density fluctuations smoothed with a given window, and let 
5g{x) be the corresponding field for galaxies of a given type. We treat them as random fields, 
both with probability densities of zero mean by definition. Denote (5^) = o"^ and (5^) = S. 
Consider the "biasing" relation between galaxies and mass to be a random process, specified 
by the conditional probability function B{6g\S). The common deterministic biasing relation, 
6g = b{S)S, is replaced by the conditional mean. 



The statistical character of the relation is expressed by the conditional moments of higher 
order about the mean, such as 



{SG\6)^b{6)6. 



(18) 



{{6a-b6)'\6)=al and {{6a - b6f\6) = S^. 



(19) 



This statistical nature of biasing leads to a different "biasing parameter" for each specific 
application. 
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Take for example the ratio of variances, = {S'^)/{S^), such as being obtained by a 
ratio of power spectra or two-point correlation functions, or by comparing the mass function 
of clusters to the variance of Sq at 8 /i~^Mpc (White, Efstathiou & Frenk 1993). One can 
prove in general that (5™) = ( (C I'^)5g and therefore, {5l) = {b'^{5) 5^) + {al{5)). Thus, 
in the simple case where b{6) is constant, 62 is an overestimate of b by 

b2 = b{l + A,Y/', A,^{al)/ia'b'). (20) 

Another common way of estimating j3 is via linear regression of the noisy field — V-t;(^) 
f{fl)6{x)] on 6a{x) (§ 1^ ), or via a regression of the corresponding velocities. The slope of 
the forward regression of 6g on 6 is bf = (ScS) / (5^), and the slope of the inverse regression of 
5 on Sg is b^^ = {55g) / {5^)- In the case where b is constant, bf = b, and bi is an overestimate, 

6, = 6(1 + A2). (21) 

The promising method of estimating f3 from large-scale redshift distortions measures yet a 
different quantity. It turns out that most methods for determining (3 lead similarly to an 
underestimate. 

The level of the effect depends on the actual values of A2 and similar parameters. One 
way to estimate the natural biasing scatter at a given smoothing scale is by investigating 
goodness of fit of the density fields of mass and light and the model of deterministic biasing. 
By requiring that = 1 P^r degree of freedom one can estimate the scatter needed in 
addition to the known errors. For example, Hudson et al. (1995) estimated for optical galaxies 
versus POTENT Mark III mass at 12 /i~^Mpc smoothing (Jb ~ 0.15, which corresponds to 
A2 ~ 0.25. Alternatively, one can estimate ah from theoretical simulations. For example, 
preliminary hydro simulations (Cen & Ostriker 1993) yield cxfe = 0.25 under 10 /i~^Mpc 
Gaussian smoothing, i.e. A2 ~ 0.4. If 6 = 1] = 1, then the (3 values derived by the various 
methods are expected to span the range 0.7 < /5 < 1, and this is solely due to the dispersion 
in the biasing relation. A large skewness in B{5g\5) may stretch this range even further. 

A relevant moral from the biasing uncertainty is that methods for measuring VL inde- 
pendent of density biasing (§ ^) are desirable. However, it has to be born in mind that 
the galaxies may also be biased tracers of the velocity field of the matter. Such a "velocity 
biasing" would affect any attempt to extract dynamical information from large-scale veloci- 
ties. The expected magnitude of the velocity biasing in the standard scenarios of structure 
formation is a matter of debate, and even it's sign is unclear {e.g., Summers et al. 1995). 
Based on recent simulations it seems likely to be limited to a ~ 10 — 20% effect. 



6. COSMOLOGICAL PARAMETERS 

The previous sections discussed measurements of the mass-density parameter VL (directly 
or via (3) from large-scale structure on scales 10 — 100 /I'^Mpc. In this section we try to put 
these estimates in a wider perspective (see Dekel, Burstein & White 1997 for a review). 
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One very interesting large-scale constraint that has not been discussed here is based 
on cluster abundance, that can be predicted for a Gaussian field via the Press-Schechter 
formalism, and is quite insensitive to the shape of the power spectrum. The current estimates 
are agfi"''^ ~ 0.5 - 0.6 (White, Efstathiou & Frenk 1993; Eke et al. 1996; Mo et al. 1996). 
This is only slightly lower than the estimates of a^^^'^ — 0.7 — 0.8 from the power spectrum 



of the peculiar velocity data (§ ^72] , § Note that this quantity is related to /3j via 

a^Vf'^ = asi/Sj, where asi is the rms fiuctuation of IRAS galaxies in a top hat window of 
8 h~^Mpc. With the estimate from the IRAS 1.2 Jy survey of ag/ = 0.69 ±0.05 (Fisher et al. 
1994a), the results from cluster abundance and from the 6 — 6 POTENT-IRAS comparison 
(§ [5. 21 ) are in pleasant agreement. 

Constraints from virialized systems such as galaxies and clusters on smaller comoving 
scales of 1 — 10/i~^Mpc {e.g., Primack 1997; Bahcall 1997; Peebles 1997) typically yield low 
values of ~ 0.2 — 0.3, but with several loopholes. Most interesting among these is the 
constraint involving the baryonic fraction in clusters from X-ray data and the estimates of Qb 
from the observed deuterium abundance and the theory of big-bang nucleosynthesis. With 
baryonic fraction in the middle of the observed range fb = (0.03 - 0.08)/i-3/2 (White et al. 
1993; White & Fabian 1995), and with the recent estimates of fife/i^ ^ 0.025 (Tytler et al. 
1996; Buries & Tytler 1996), the current estimate is ~ 0.5/ig5 . This result favors a low 
value of fl, but = 1 cannot be definitively excluded. 

The global cosmological measures commonly involve combinations of the cosmological 
parameters, such as the mass Q (= Qm) and the contribution of the vacuum energy density 



^A- Constraints in the Q — Q\ plane are displayed in Figure 15, and briefiy discussed below. 



Occam's Razor. The above working hypotheses, and the order by which more specific 
models should be considered against observations, are guided by the principle of Occam's 
Razor, i.e., by simplicity and robustness to initial conditions. It is commonly assumed that 
the simplest model is the Einstein-deSitter model, = 1 and Q\ = 0. One property that 
makes it robust is the fact that Q remains constant at all times with no need for fine tuning at 
the initial conditions. The most natural extension according to the generic model of infiation 
is a fiat universe, fltot = 1, where Q can be smaller than unity but only at the expense of a 
nonzero cosmological constant. 

Classical Tests of Geometry. The parameter-dependent large-scale geometry of space- 
time is refiected in the volume-redshift relation. There are two classical versions of the tests 
that utilize this dependence: magnitude versus redshift (or "Hubble diagram") and number 
density versus redshift. The luminosity distance and the angular-diameter distance to a 
redshift z, which enter these tests, depend on Q and Q\. At 2; ~ 0.4, these distances happen 
to be (to a good approximation) a function of the combination Q — Q\ (not go) (Perlmutter 
et al. 1996). 

The main advantage of such tests is that they are direct measures of global geometry. 
Supernovae type la are the popular current candidate for a standard candle, based on the 
assumption that stellar processes are not likely to vary much in time. 

The first 7 supernovae analyzed by Perlmutter et al. (1996) at 2 ~ 0.4 yield —0.3 < 
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f2 — f^A < 2.5 as the 90% two-parameter likelihood contour (Fig. 15). For a flat universe 
they find for each parameter Q = 0.941q28, and < 0.51 (or equivalently Q > 0.49) at 
95% confidence. 

Number Count of Quasar Lensing. This is a promising new version of the classical 
number density test. When ^a is positive and comparable to Q, the universe should have 
gone through a phase of slower expansion in the recent cosmological past, which should be 
observed as an accumulation of objects at a specific redshift of order unity. In particular, 
it should be refiected in the observed rate of lensing of high-redshift quasars by foreground 
galaxies (Fukugita et al. 1990). The contours of constant lensing probability in the VL — ^a 
plane for ~ 2 happen to almost coincide with the lines VL — f^A = const. The limits from 
lensing are thus similar in nature to the limits from SNe la. 

This test shares all the advantages of direct geometrical measures. The high redshifts 
involved bring about a unique sensitivity to ^Aj compared to the negligible effect that ^a 
has on the structure observed at z <^ 1. 

From the failure to detect the accumulation of lenses, the current limit for a flat model 
is ^A < 0.66 (or n > 0.36) at 95% confidence (Kochanek 1996) (Fig. |15D. 

Microwave Background Acoustic Peaks. This test is expected to provide the most 
stringent constraints on the cosmological parameters within a decade. The next generation of 
CMB satellites (MAP, to be launched by NASA in 2001, and in particular Plank, scheduled 
by ESA for 2004) are planned to obtain a precision at ~ 10 arc-minute resolution that 
will either rule out the current framework of GIfor structure formation or will measure the 
cosmological parameters to high precision. Detailed evaluation of Plank shows that nominal 
performance and expected foreground subtraction noise will allow parameter estimation with 
the following accuracy (ignoring systematics): Hq ± 1%, Qtot ± 0.005, ^a ± 0.02, Qf, ± 2%. 

Current ground-based and balloon-born experiments provide preliminary constraints on 
the location of the first acoustic peak on sub-degree scales in the angular power spectrum of 
CMB temperature fiuctuations, /(/ + l)Ci. In the vicinity of a fiat model, the first peak is 
predicted at approximately the multipole Ipeak — 220(fi + ^a)^^^^- The results of COBE's 
DMR (/ ~ 10) provide an upper bound of fi+i^A < 1-5 at the 95% confidence level for a scale- 
invariant initial spectrum (and the constraint becomes tighter for any "redder" spectrum, 
n < 1) (White & Scott 1996). Several balloon experiments (/ ~ 50 — 200) strengthen this 
upper bound {e.g., the Saskatoon experiment, Scott et al. 1996). The Saskatoon experiment 
and the CAT experiment (/ ~ 350 — 700) yield a preliminary lower bound of f2 + ^a > 0.3 
(Hancock et al. 1996) (Fig. |15|). 

The Age of the Universe. Measured independent lower bounds on the Hubble constant 
and on the age of the oldest globular clusters provide a lower bound on -ffo^o (= 1.05ht, 
where Hq = 100/ikms~^Mpc~^ and to = lOtGyr), and thus an interesting constraint in the 
fl — fl\ plane. The exact expressions are computable in the various regions of parameter 
space. A useful crude approximation near -ffo^o ~ 2/3 is — O.T^a — 5.8(1 — l.Sht). 

Progress has been made in measuring Hq via the HST key project detecting Cepheids 
in nearby clusters for calibration of TF distances, and via accurate distances to SNe Type 
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Fig. 15. — Current limits (~ 2a) on the cosmological parameters and Q\ from global measures: 
luminosity distance of SNIa, lens count, the location of the CMB peak, and the age versus Hubble 
constant. The short marks are the one-parameter 95% limits from SNIa and lenses for a flat 
universe. Also shown (vertical line) is the 95% lower bound on 0, from cosmic flows. The most 
likely value of lies in the range 0.5 to 1. The Einstein-deSitter model is permitted. An open 
model with ~ 0.2 and Q\ = 0, or a flat model with 0, ~ 0.3 and J7a — 0.7, are ruled out (Dekel, 
Burstein & White 1997). 

la. The new calibration of local Cepheids by the Hipparcos astrometric satellite (Feast & 
Catchpole 1997) seem to have reduced the estimates of Hq by ~ 10%. The indications from 
SNe velocities for a local void of radius ~ 75 h~^Mpc out to the Great Wall (Zehavi et al. 
1997), takes another ~ 5% from Hq as measured by TF distances from within the void, and 
brings the various estimates into agreement at/i~0.6±0.1. 

The Hipparcos calibration of the distances to local subdwarf stars (Reid 1997) had an 
even more dramatic effect on the estimates of the ages of the oldest globular clusters {e.g., 
Van den Berg et al. 1996). The current estimates seem to be t ~ 1.2 ± 0.2 {e.g., M. Bolte, 
private communication). Thus, the most likely value of HqIq is not far from 2/3, consistent 
with the standard Einstein deSitter model. 



6.1. Conclusion 

We conclude with a summary of the main implications of the observed cosmic flows. 

Gravitational Instability. The strongest evidence for gravitational origin of structure 
comes from the growth rate of fluctuations as indicated by the comparison of the 5T/T~ 10~^ 
fluctuations at the last scattering surface and the ~300kms~^ motions over ~ 100 /i~^Mpc 
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scales in our local neighborhood. 

Initial fluctuations and Dark Matter. The COBE measurements of CMB fluctuations 
at large angular scales and the comparison to the observed flows indicate a power spectrum 
near scale invariance, n~l. 

The bulk velocity in a sphere of radius 50 h^^Mpc about the LG is V50 = 375±85 kms~"^. 
The mass power spectrum deduced from the peculiar velocities has an amplitude of Pq.i^^'^ — 
(5 ± 2) X 10^{h-^Mpcf at 0.1{h-^Mpc)-^. This extrapolates to crgll"-^ = 0.8 ± 0.2 on 
smaller scales. 

For COBE-normalized CDM models, a likehhood analysis of the mass power spectrum 

yields Vtn^hQ^ = 0.7 ±0.2. A comparison to preliminary detections of the first acoustic peak 
in the CMB angular power spectrum requires that n >, 0.9, and that Qf, ~ 0.1. Thus, within 
the family of CDM models, most successful in matching the current LSS data are either of 
the following variants: (a) fl — 1 with a tilted spectrum n ~ 0.9, (b) ft ~ 0.5, with or 
without a cosmological constant, and n = 1, and (c) Q = 1 with 20% hot dark matter. A 
high baryonic content (and relatively low Q) may be required to explain a peak in the galaxy 
density power spectrum at ~ 125/i~^Mpc, if confirmed (Broadhurst et al. 1990; S. Landy 
based on LCRS, private communication; Cohen et al. 1996; Einasto 1997). 

In view of the nucleosynthesis constraints on the baryonic density, the high Q indicated 
by the motions requires non-baryonic dark matter. The observed mass-density power spec- 
trum on scales 10—100 /i~^Mpc does not yet allow a clear distinction between the competing 
models involving baryonic, cold and hot dark matter and possibly a cosmological constant. 
I do not think that any of the front-runner models can be significantly ruled out based 
on current tests, contrary to occasional premature statements in the literature about the 
"death" of a certain model. I predict that were the dark matter constituent (s) to be securely 
detected in the lab, the corresponding scenario of LSS would find a way to overcome the 
~2(T obstacles it may be facing now. 

Galaxy Biasing and j3. Generally speaking, galaxies trace mass. For each of the different 
smoothing scales, the data of velocities and redshift surveys are consistent with Gl and linear 
biasing (properly modified in the tails). However, the best estimates of /3j span the range 
0.5 — 1.0. This can be explained by the fact that, when inspected in detail, the biasing 
scheme involves scale dependence, non-linear features, and intrinsic scatter. It is difficult to 
distinguish non-linear biasing from non-linear gravitational effects. 

The Value of Vt. Methods based on virialized objects tend to favor low values of f2 ~ 0.2, 
but with plausible loopholes. 

The current peculiar- velocity data provide in several different ways a significant (> 2a) 
lower bound of > 0.3. This bound is independent of A, Hq, and the biasing relation between 
galaxies and mass. The range of (3 values obtained on different scales by different methods 
may be partly due to underestimated errors and partly due to non-trivial biasing. 

The global measures of geometry provide a lower bound of similar nature, ^2 — > 0.3. 
The age constraints, which used to favor low values of Vt until recently, seem to agree with 
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f2 ~ 1 according to the new calibration of the distance scale by Hipparcos. 

The data is thus consistent with Q = l. Based on the whole range of constraints, and 
ignoring the Occam's razor desire for simplicity, the most likely value may be argued to be 
Q ~ 0.5. Values of Q = 0.3 and below are significantly ruled out. The data are consistent 
with the general predictions of Inflation: fiat geometry and Gaussian, almost scale-invariant 
initial fluctuations. 
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